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Abstract. Quantum entanglement plays a crucial role in quantum information, 
quantum teleportation and quantum computation. The information about the 
entanglement content between subsystems of the composite system is encoded in the 
Schmidt eigenvalues. We derive here closed expressions for the spectral density of 
Schmidt eigenvalues for all three invariant classes of random matrix ensembles. We 
also obtain exact results for average von Neumann entropy. We find that maximum 
average entanglement is achieved if the system belongs to the symplectic invariant 
class. 
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1. Introduction 

Quantum entanglement serves as the measure of quantum non-local correlations between 
the subsystems of a composite system. It has attracted a great deal of attention in recent 
times because of its crucial role in quantum information theory, quantum computation 
and quantum teleportation [U [2]. Quantum entanglement has bizarre consequences 
which have baffled physicists over the years [3llll[5]. Despite this, quantum entanglement 
is now universally acknowledged as a useful resource which can be harnessed for a number 
of apphcations, e.g., superdense coding [6l [7], quantum teleportation [8], quantum 
cryptography [9], creation of a quantum computer [10] etc. 

In most cases it is desirable to have maximum possible entanglement between the 
subsystems of the compound system. It is known that random pure states lead to 
nearly maximal average entanglement. Most of the work so far has concentrated on 
random pure states which belong to the unitary invariant class (/3 = 2) jTTl [131 
ttH [ISl UHl HTJ HE], which corresponds to systems with broken time- reversal symmetry. 
It is natural to investigate time-reversal symmetric systems, i.e., random pure states 
belonging to the orthogonal (/3 = 1) or symplectic (/3 = 4) invariant classes. Recently 
there has been some work on systems with orthogonal invariance and now there are 
some important results available for systems with finite as well as large Hilbert space 
dimensions [I9l [201 EH [22]. For systems with symplectic invariance there has been 
work concerning only large dimensions and for extreme eigenvalue statistics [191 1231 [23] . 
For these systems the distribution of Schmidt eigenvalues and various entanglement 
measures have not been investigated for finite dimesional cases. 

Our purpose in this paper is to treat the random pure states belonging to all three 
invariant classes of random matrices [251 EE] on equal footing and obtain results valid for 
systems with arbitrary Hilbert space dimensions. We derive here the spectral density of 
Schmidt eigenvalues for arbitrary dimensions. We also derive the average von Neumann 
entropy and show that maximum average entanglement is achieved if one considers 
systems belonging to the symplectic invariant class. 

The presentation scheme in this paper is as follows. We begin with the brief 
introduction to the Wishart-Laguerre ensemble and give exact results for the spectral 
density in Section 2. In section 3 we briefly discuss the quantum entanglement problem 
in bipartite systems. Section 4 deals with the flxed-trace Wishart-Laguerre ensemble 
which is relevant for the entanglement problem. We obtain in this section closed results 
for the spectral density of Schmidt eigenvalues for all three invariant classes of random 
matrix ensembles. In section 5 we calculate the average von Neumann entropy. By 
comparing the von-Neumann entropies in the three cases, we show that, as far as average 
value is concerned, maximum entanglement is achieved for /3 = 4. Finally, we conclude 
in section 6 with summary and some general remarks. Some proofs of the results are 
outlined in the appendices. 
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2. Wishart-Laguerre ensemble 

Wishart-Laguerre ensembles were introduced by Wishart [27] in connection with the 
analysis of multivariate distributions. These ensembles comprise the A?^- dimensional 
matrices H = XX^ where the X are N x M dimensional matrices having its elements 
as independent and identical (iid) Gaussian random variables. 

Wishart-Laguerre ensembles arise explicitly in a number of problems in completely 
unrelated areas. Besides the quantum entanglement problem [2HI EHl [30] some other 
examples are chaotic mesoscopic systems [31], low energy QCD and gauge theories [32] , 
quantum gravity [33l[3l], quantitative finance [35], [36], communication theory [371 [38]. 
gene expression data analysis [391 HO] etc. 



2.1. Joint probability density of eigenvalues 

Consider the elements of the matrices X constituting the Wishart-Laguerre ensemble 
XX^ as real, complex or quaternion-real zero-mean iid Gaussian variables. These three 
cases lead to the three invariant classes of ensembles referred to as orthogonal, unitary 
and symplectic ensembles and are designated by the Dyson index /3 = 1, 2 and 4 [25|[26]. 
We choose N < M for definiteness. The probability distribution followed by X is 

P(X)ocexp(-^), (1) 

where is the variance of each distinct real component of matrix elements of X. The 
respective joint probability density (jpd) of eigenvalues {xj G [0,oo),j = 1,...,A^) of 
XX^ for the three cases come out as [261 HD 112] 

N N 

P'^'\{^}) = Cg^ |A^({x})|%xp ( - , (2) 

i=i j=i 

where {x} denotes {xi, ...,xn}, ^n{{x}) = Y[j<k(-'^j ~ ^k) is the Vandermonde 
determinant and the parameter uj is defined as 

uj = ^{M -N (3) 

The normalization C^^^ ^ can be evaluated using the Selberg integral |i25j and turns out 
to be 

^M,N llr(|, + i)r(|(,_i) + ^ + i)- 

Note that ([2]) is obtained from ([T|) for the choice v"^ = 1/2. Some authors choose f ^ = 1 
(see for example pS]), in which case apart from the change in normalization, the jpd in 
([2|) will have additional factor of 1/2 in the exponential. 
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2.2. Spectral densities for eigenvalues 

In many cases one is interested in the study of quantities which are hnear statistics on 
the eigenvalues. These quantities do not contain products of different eigenvalues. The 
averages of these can be evaluated using a single integral over the level-density of the 
eigenvalues, defined as 

dx2-..dxN. (5) 



From random matrix theory [25j we know that, for the classical ensembles, the above 
density can be obtained in terms of classical orthogonal polynomials for /3 = 2 and in 
terms of classical skew-orthogonal polynomials for /3 = 1, 4 [HI [45], |46l HTJ l48] . 

For the Wishart-Laguerre ensemble defined by jpd ([2]), the exact expressions for 
level-densities are known for all M, [49j in terms of associated Laguerre polynomials 
L^ll\x) [50] and the corresponding weight function Wh{x) = x^e~^, b > —1. Consider 
the parameter a defined by 

2a + l = \M-N\. (6) 

We have for unitary ensemble (/3 = 2), 



N-l 



= -...i(x)^;^^^^^ [L;j-^)(a:)Lr^)(x) -Lr^)(x)Lr_f^(a:)] .(7) 

The second line in the above equation has been obtained using the Christoffel— Darboux 
sum [So]. Other expressions equivalent to the above sum are also possible; see (IC.ll) . 

For orthogonal ensemble (/3 = 1) the level density as well as higher order correlation 
functions are given in terms of skew-orthogonal polynomials These skew-orthogonal 
polynomials can be expanded in terms of orthogonal polynomials jlTJ HHj. Closed form 
expressions for the level density can be obtained using these expansions after some 
tedious algebra. We obtain for even N [19], 

<Ln(A^;«;a;) = 2i?p)(iV;a;2x) 

iV-2 

1^ w(2a+l)/D \ r(^— )r(/i+l) (2a+l)/^ n 

-^,.,,(2x)4_, \2x) g r(/. -T^'^^-^^ ^'"^ 

, , xW2a+i).^ X r(^) r2r(a + l,x) 1 , . 

+ti;,(x)LV_i (2a;)^7|Mpiy [^X^TTT " J' 
while for odd we get 
<odd(^;«;a;) = 2i?f)(iV;a;2x) 

N-l 

^w.{x)Lf^^\2x)^^^L_^ (9) 
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The above two expressions may be put together as a single expression by defining 

c = A^(mod2), (10) 
so that we have for all A^, 
R^^\N;a;x) = 2R^^\N; a;2x) 

N-2+C 

Nr(2a+l)/„ N r ) r(/i + 1 — |) r(2a+l) /r, x 

-w,.,,{2x)L,_, (2x) g 22.r [-±f±l) r{f^ + a Ti^^^^^^-^'"^ 
+w.{x)L^^^^\'\2x)^^^^[{l - c)^^!^^ + (2c - 1)]. (11) 

Here T{b,x) represents the incomplete Gamma function, T{b,x) = y^~^e~ydy. 

Finally we consider symplectic ensemble (/3 = 4). Similar to /3 = 1, the level density 
in this case also is given in terms of the corresponding skew-orthogonal polynomials. We 
use the expansion of these skew-orthogonal polynomials given in terms of orthogonal 
polynomials to obtain the following closed expression for the level-density [15] : 

R^^\N; a; x) = ii?S^^(2A^; 2a + 1/2; x) 

I xr(4a+2)/ nV^^ r(A^ + l)r(/i + i) r(4a+2)/ x 

4a+2v ; 27V V ^ 24«+3r(A^ + 2a + |)r(/i + 2a + 2) v ; v ; 

Note that in view of definition given by ([6]), the parameter a assumes only integral 
or half-integral values (the Wishart case), however, the above results hold for all real 
a > — 1 (the Laguerre case). 

Comparison between the above theoretical results and corresponding numerical 
simulations is shown in figure 1. They are in excellent agreement. 




2 4 6 



10 12 14 16 18 



Figure 1. Density, p('^)(iV; a; x) = R^f\N;a]x)/N, of eigenvalues for /3 = 1,2,4. 
Solid lines are the theoretical predictions, whereas the symbols represent the simulation 
results. 
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For large A^, the density is described by the Marcenko-Pastur formula [5T| [52] , 



0, otherwise. 



i?f ^iV; a; x) = <^ 7r/3x ' " " +' (13) 



Here x± = (/3M/2)(1 ± ^N/Mf. 
3. Bipartite system 

Consider a bipartite partition of an iVM- dimensional Hilbert space 'H'^^^^^ consisting 
of subsystems A and B which span Hilbert spaces •Hf ) and -Hf ) such that -H^^*^) = 
Hf^®n'i^\ For example, A may be a set of spins and B can be a heat bath. Without 
loss of generality we choose M > N . We take \^) to be a quantum state of the composite 
system and \(y^) as the complete basis of "H^^^ and "Hg^^ respectively. The state 
|\E') can then be expanded as the linear combination 

N M 
i=l a=l 

We focus on the situation where |\E') is a pure state, so that the density matrix of 
the composite system is given by p = |\E')(\E'|, with the constraint Tr[p] = 1, which is 
equivalent to (^fl^f) = 1. 

The density matrix for the pure state \^) can be written as 

N M N M 

p = ® (15) 

i=l a=l j=l 7=1 

The reduced density matrix for subsystem, say A, can be obtained by tracing out the 
subsystem B. We obtain 

N N 

Pa = EE^mK^)(/I, (16) 
i=i j=i 

where Wij = J2a=i di^aQj^a can be viewed as the matrix elements of some matrix 
W = GG\ G being N x M dimensional. In the diagonal basis, ([T6|l can be written as 

N 

PA = 5^A.|Af)(Af|. (17) 

1=1 

The Aj are referred to as Schmidt eigenvalues. The trace condition restricts the Schmidt 
eigenvalues to the following constraint 

N 

^A, = tr(GG^) = 1. (18) 

i=l 

|\E') is called a random pure state if the coefficients Qi^a are chosen at random from some 
given probability distribution. 
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4. Fixed trace Wishart-Laguerre ensembles 

As discussed in the previous section, the Schmidt eigenvalues correspond to Wishart type 
of matrices. However, the conservation of probabihties for quantum states constrains 
the trace of these matrices. In earher works Wishart-Laguerre ensembles have been used 
to model random quantum states [29l |30]. However, they do not provide a complete 
quantitative description of the random quantum states due to the neglect of fixed trace 
condition [30]. The deviations are prominent particularly when the Hilbert space- 
dimension of the smaller subsystem is not large. Detailed discussion of superiority of 
fixed-trace ensemble over the ordinary Wishart-Laguerre ensemble for describing random 
quantum states can be found in [T71 [18] . 



4.I. Joint probability density of Schmidt eigenvalues 

Consider the matrix elements Qi^a as iid (real, complex or quaternion-real) Gaussian 
variables with zero mean, such that the fixed trace condition f|T8|) is satisfied. The 
probability distribution for G is given by 

V{G) oc 6{tT{GG^) - 1). (19) 

Correspondingly the jpd of Schmidt eigenvalues (Aj G [0, = 1, A^) is obtained as 

N N 

^^^^{A}) = 0(E^^- \^Nm)fi[^i (20) 

i=l 0=1 

Note that, with the fixed trace condition, (|T8l) . the exponential term in (jj]) becomes a 
constant and one is left with f l20|) . The normalization constant in the above equation is 

where C^^^ is given by (jl]); see Appendix A The three /3 values (1, 2, 4) in (j20ll 
again correspond to the cases when the coefficients gi^a in dlSD are real, complex or 
quaternion-real. Also, the choice of variance of Gaussian elements is immaterial in this 
case. 



4-2. Spectral densities for Schmidt eigenvalues 
The level-density of Schmidt eigenvalues is, 

/•oo /"OO 

nf\M,N-\) = N ■•■/ p(^)(A,A2,...,Aiv)c/A2...rfAjv. (22) 
Jo Jo 

Note that because of the presence of delta function constraint the original limits of 
integration [0, 1] can be replaced by [0, oo). 

We give here the exact expressions for the level density of Schmidt eigenvalues 
obtained from the jpd (l20l) . Proof is briefly outlined in Appendices B and C. 
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We need the following to express the desired results compactly: 
, (-l)^r(|m + l)r(|mn) 

^ ~ |r(/i + l)r(|n - /i)r(|(mn - m + n) - /i - l)r(/i + |(m - n) + 2) ' 

Tjm,n ^ I ^) ^\2)'-y2)^y2) (24) 

^ vri/2r(H^)r(/i + m-n + l)r(/^ + l)r(n-/i)' ^ ^ 



om— np / m+]_ \ p / mn N 

p)m,n ^ I 2 V 2 ^ /of^N 



2 V 2 ' 



22/.-2iV/^2m-2n + l)r(2n-z/)r(/i + i)r(u + m-n+i) , , , 



r(m + i)r(n + i)r(2/i - z/ + 1) 



X (r^ - ^n, /i - |(m + l)(n - 1); 1 + |(m - n); , (28) 

m, n; y) = AJ'"'" [|n ^,1,1; m, n, - (^^n - - l) 0, 1; m, n; y)] , (29) 

W^(m,n;y) = (2y) ^+"^-"+1(1 - 2y)-'^+"i"— +""1 

^ / 2 TflTl 1 \ 

X 2^1 (^l,--(m-n- 1); — - m + n-/^; 2- - j. (30) 

Here 2-F'i(^, "O^; y) = 2-^1 (^ "O^; n; y)/T{c) is the regularized hypergeometric function. 

We first consider the unitary case (/3 = 2), which has been considered by several 
authors. In [T7j the result was announced for the square case (M = A^) in terms of sum 
involving hypergeometric 3F2. Later the result was extended to rectangular case and was 
given in terms of sum comprising hypergeometric 5F4 [18]. In [21] a triple sum expression 
has been given. We give here a result, valid for all M > N, which is much simpler than 
those obtained earlier and is given as a single sum over the hypergeometric 2-^1- We 
have 

N-l 

nf^ (M, N;X) = Y, %i2; M, N; A), 

11=0 

(31) 

The spectral density for the orthogonal case (/3 = 1) has been obtained for N even 
in [21]. We consider here both even and odd cases. For even N, we obtain 

^Sven(M, AT; A) = [4 5^ r, (4; -, -; 2a) + Bff'^'U.iM, N, A) 

11=0 ij.=0 



N-l 



— D^^'^S~M+N~i (a, 1, 2; ^r', A 



' ' ' 2 ' 2 ' " 



(32) 



while for odd we have, 



Af-l 



<L(M,iV;A)= [4 5^r,(4;^,^;2A 



fi=0 

JV-l 
Af— 1 ~2~ 

-EEC^.r5.(4,iV-2.,l;f.f;2A 



e 



(33) 

Similar to flTT]) the two expressions for even and odd cases may be combined into one 

as, 

N-l N-l 

TZ? (M, N- A) = [4 5^ r, (4; - , -; 2A) + (1 - c) ^ A) 



JV-l 



M=0 

N-2+C 



/x=0 



E E C^,r'5.(4,iV-2.-l + c,l;f,f;2A 

^=0 u=0 



el--A 



-l)^^£)^"'^^>S -A/+jv-i ( 4, 1, 2; — , A 



(34) 



Q{y) in the above equations is the Heaviside theta function, being 0, 1 respectively for 
y < and y > 0. 

Finally, we consider the symplectic case (/3 = 4). Note that we do not take into 
account the Kramers degeneracy explicitly and use the jpd given by fl2Ul) . We obtain 

N-l 2)1 

n['^ (M, N;\)=Y, T(4; M,N-\)-Y,Yl ^mT'5.(4, 0, 1; M, N; A). (35) 

/x=0 fj.=0 u=0 

Of particular interest is the case of = 2 and arbitrary M > 2, as it corresponds 
to a physical situation where a single qubit is entangled to a heat bath. In this case the 
expression for level density is simple and can be obtained directly from (|20|) . We have 
for all three P, 

nf\M, 2; A) = «:(/3'A'^)Af (^^-i)-i(l - X)-^^''-^^-^\l - 2A|^ 

where k*^^'*^) is given by 

''2*^-2(M-l), 13 = 1, 

(M-l)r(2M) 



(36) 



K 



W,M) 



(r(M))2 ' 

(2M- l)(2M-2)r(4M) 
6(r(2M))2 ' 



(3 = 2, 
(3 = A. 



(37) 
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For large M, A close to 1/2 dominates the level density (1361) and we obtain, 

n<f\M.2;X, « _ 2A|^. (38) 

As can be seen in figure [2] as M increases the two maxima move closer and closer towards 
A = 1/2 and become sharply peaked like delta function for M oo. 

We have numerically generated fixed trace Wishart-Laguerre ensemble by 
considering matrices XX^ /ti(X X^), where XX^ as above form the Wishart-Laguerre 
ensemble. Comparison of the above theoretical results with numerics is shown in figure. 
They agree perfectly. 




@H 0.25 0.5 0.75 0.25 0.5 0.75 0.25 0.5 0.75 1 




0.25 0.5 0.75 0.2 0.4 0.6 0.8 01 0.25 0.5 0.75 1 



Figure 2. Density, p[f^(M,7V;A) = n[^\M, N; X)/N , of Sclimidt eigenvalues for 
/3 = 1,2,4 for several M, N values. Solid lines are the theoretical predictions, whereas 
the symbols represent the simulation results. 

The level densities obtained above give good description of the Schmidt eigenvalues. 
These expressions are important because the statistical properties of entanglement are 
encoded in the Schmidt eigenvalues. For the unitary case the one-level density has been 
found to describe accurately the distribution of the Schmidt eigenvalues over the whole 
range for coupled quantum systems exhibiting chaos |T7j. This confirms the existence of 
universality for the distribution of the Schmidt eigenvalues in quantum chaos. The level- 
densities can be used to calculate the averages of quantities which are linear statistics 
on Schmidt eigenvalues, e.g., purity, von Neumann entropy, Renyi entropy etc., thereby 
giving insight into the extent of entanglement. 
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5. Average von Neumann Entropy 

As discussed above, the information about the degree of entanglement is stored in 
the Schmidt eigenvalues. One can construct a suitable measure of entanglement by 
considering a function of these eigenvalues whose value will tell us how entangled a 
pure state is. Von Neumann entropy serves as an important measure of degree of 
entanglement. It is an extension of classical entropy concepts to the field of quantum 
mechanics. In terms of the Schmidt eigenvalues, it is defined as 

N 

= -^A^lnAj-. (39) 
i=i 

It acquires the maximum value of In A, when each of the Scmidt eigenvalues assume 
the value 1/A. On the other hand the minimum value of is obtained when one of 
the eigenvalues is one and rest are zero. These two scenarios correspond respectively to 
maximally entangled and completely unentangled (separable) cases. 

To calculate the average von Neumann entropy, E^^^ = {£^^^), we need to perform 
the following integral: 

poo poo ^ 

E^^\M,N) = - / ■■■ / ( VAanA,)p(^)({A})dAi...rfA^. (40) 
^0 Jo ^ 

As in (l22l) . we have again considered the limits of integrations as [0, oo) instead of [0, 1]. 
Von Neumann entropy being a linear statistic, the symmetry of the eigenvalues in jpd 
allows us to reduce the above average involving A integrals to an average involving a 
single integral, i.e., 

POO 

E^^\M,N) = - WnX n^^\M,N;X)dX. (41) 



^0 

The above calculation can be done using the results for level density given in preceding 
section. However, it turns out that the above average over the fixed trace ensemble 
can be reduced directly to an average over the ordinary Wishart-Laguerre ensemble by 
introducing an auxiliary gamma function integral [12]. The latter route to calculation 
of average von-Neumann entropy is somewhat more straightforward and we follow it 
below. 



As shown in Appendix D we have 



E(« (M, N)=i; + l) - ^ l^^x In x R^^ (A; a; x) dx. (42) 

where R^^"^ (A; a; x) is the level density for ordinary Wishart-Laguerre ensemble, as given 
in ([5]) and tplx) represents the digamma function defined by 

1 f°° 

^(y)= e-W^-Hnrdr. (43) 

r(2/) Jo 
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Equation fl42p gives a compact expression for the average von Neumann entropy for all 
three P values. For M > >> 1, the average von Neumann entropy can be found from 
I2l) using the Marcenko-Pastur density given by ( |T3l) . Also in this limit 



(44) 



Using ([I3D, the integral in (02]) is obtained as {(3MN/2)\n{(3M/2) + (3N'^/4:, thereby 
giving the large result for all three /3 as 

(45) 



E'^^\M,N) = \n{N) - 



2M 



The above result was obtained for /3 = 2 in [12]. See [29] for /3 = 1. 




Figure 3. Average von Neumann entropy for (a) N = 2,3,4 and M = 2 to 16 (b) 
M = 8 and 2 to 8. 



Exact closed form expressions for the average von Neumann entropy for arbitrary 
M, A^ can be calculated using fH2|) and the results in section l2T2l To evaluate the integral 
in fH2|) we use the following identity: 



/CO g i>oo 

x^\n{x)f{b,x)dx = — J x^f{b,x)dx 



x'^J^dx. 



db 



(46) 



The evaluation of the above integral for the three invariant classes involves the use of 
properties followed by Laguerre polynomials [50] and its expansion (IC.2p . 

The average von Neumann entropy in /3 = 2 case is already known [T2l [T3] . We 
have, 

E<^\M.N)= l-^. (47) 

For /3 = 1, even A^, we obtain 

'MN 



2M 



, , (N-1) (M-i)(Ar-i), 2^-ir(^)r(^^ 



2M 



MN 
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iV-2 ^ 

/i=0 i/=0 



+ M-N+ l)r(/x + l)7/>(z/ + M - AT + !)(-;/ - 1)2^+1 

r(/x + ^^i^)r(iv - u)T{iy + i)r(2/x + 2) 



X 



2'^r(^)r(^) 

M-N+d 

"^"^ (-l)"r(/i + I. + M^i^t2^)(V'(/i + I. + M^i^t2^) _ In 2) 



,=0 2Mr(z. + i)r(iv-z/)r(z. + M-iv + i)r(/i + i^) 



^^v^ ^ r(iv-z/)r(z/ + M-A^ + i)r(z/ + i) 



2^^r(^)r(^) (-2rT(. + M^) 

^ ^ Mivy? f^^r{N-u)r{u + M-N + i)r{u + i)' ^ ' 

Here, 

= ^-^^ (49) 
represents the Pochhammer function, 

rf = mod(M- A^,2), (50) 

and 

T(fc) 



d ( r(2a + 2) ^ , 3 



da V22"r(a + 2) , , — > ; • (51) 

For /3 = 1, odd A^, the result is comparatively simple. We have 

, (MN \ i2M + 2N - MN -2) , , (M\ 
E^^\M,N) =ij( + 1) + ^ '-ij(M)-^ -i^i — ] 

(M-l)(iV-l) _ (iV^-iV + 4) 
MN ^ ' 2MN 

(Am) r (im) ^1 (_2)'-r [v + ^(^ + M - iV + 1) 

Finally, for /3 = 4 we have 

E(4)(M,iV)= ~~ ^^^M^^ +2^^-^r(M)r(iv) 

V (~^)'^('^ + 2M - 2iV + - + 2M - 2iV + 2) 

^p^^ 22Mr(/i + i)r(z/ + i)r(2iv + i-z/)r(/i + M-Ar + i) 

The above expressions can be easily implemented in symbolic manipulation software 
packages like Mathematica. 

For large M, iV, ^ derives from (HTj), (iHl), (ESI) and ([53]) in the following way. 
In (jlTj) the summation term gives InA^ thereby reproducing fHSj) . In f HHj) the first four 
terms along with (jHl) give back (H5i) . In (152|) the first five terms give ( l45l) in the above 
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71 /T 

M 


A T 

i\ 


O 1 


(3=2 


(3 = 4: 


2 


2 


0.27718 


0.33333 


0.37063 


3 


2 


0.39772 


0.45000 


0.48099 


3 


3 


0.59987 


0.66230 


0.69918 


4 


2 


0.46451 


0.50952 


0.53490 


4 


3 


0.71324 


0.76988 


0.80167 


4 


4 


0.86452 


0.92240 


0.95494 


5 


2 


0.50680 


0.54563 


0.56692 


5 


3 


0.78503 


0.83490 


0.86210 


5 


4 


0.96152 


1.01441 


1.04329 


5 


5 


1.08065 


1.13262 


1.16108 


6 


2 


0.53595 


0.56988 


0.58815 


6 


3 


0.83442 


0.87844 


0.90204 


6 


4 


1.02846 


1.07596 


1.10143 


6 


5 


1.16366 


1.21165 


1.23744 


6 


6 


1.26130 


1.30789 


1.33298 



Table 1. Average von Neumann entropy evaluated to 5 decimal places for several 
M > N values. 

limit. Similarly in fl53p the first two terms lead to f l45p . For (3 = 1 and 4 the remaining 
summation terms become negligible in this limit. 

The average von Neumann entropy for the three invariant cases is shown in figure [31 
Table I displays the numerical values for several M, N values. It is clear that maximum 
average entanglement is obtained for (3 = 4. 

If for P = 4 we consider the Kramers degeneracy explicitly, the result for average 
von Neumann entropy will change to 

E^^)* (2M, 2 A^) = E^^^ (M, N)+\n2. (54) 

In this case the Hilbert spaces Ha and Hb are of dimensions 2N and 2M respectively 
with sum of all 2N eigenvalues being 1. Again we find E^^>{2M, 2N) > E^^'>{2M, 2N) > 
EW(2M,2Ar) and E(^)*(2M, 2A^)/2 > E'^^\M,N) > E^^'>{M,N). 
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6. Conclusion 

We have given exact results for spectral density and average von Neumann entropy for 
the random pure states belonging to one of the three invariant classes of random matrix 
ensembles. We have shown using von Neumann entropy measure that, as far as the 
average is concerned, maximum entanglement is achieved if one generates states from 
the symplectic invariant class. Several ways of generating random pure states have been 
proposed by various authors [29l |30l |53l EH ESl [561 EZ] • It will be of interest to test the 
theoretical predictions presented here using these and newer methods. 
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Appendix A. Relation between normalization constants in 04]) and ( 12 Ih 



First we find out the ratio between the normalizations in the jpd's given by ([2]) and 
(|20|). We have from (|20|). 



We introduce an auxiliary variable t in the above expression and define G{t) as. 



so that {Cjl^j j^) ^ = G{1). Taking the Laplace transform {t — )■ s) of G(t) we get 




(A.l) 




(A.2) 




(A.3) 



We now substitute sXj = Xj and obtain, 



-1 



(A.4) 



using fl2T]) . Using the inverse Laplace transform result (s — )■ t) 




t 



(A.5) 




we obtain 



t 



MN-l 




2 



(A.6) 



Finally, for t = 1 we get 



V 2 A 
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(A.7) 



Equation (lA.Tp may also be obtained using the auxialiary gamma function integral 



method discussed in [Appendix D 



Appendix B. Relation between level densities in ([5]) and l \22\} 

The level density for the Schmidt eigenvalues is (we drop the arguments M, N in the 
following) 

/■oo poo ^ ^ 

n[^\\^) = N / dX2--- rfA^C^5( VA,— l)|A^({A})|mA^ (B.l) 

"^0 -^0 S=l ^ k=l 

We introduce the auxiliary variable t and consider R^'^'*(Ai; r) defined by 

poo poo ^ ^ 

Rf\\r,t)=N d\,--- rfA;vC^5( VA,-t)|A^({A})|mA^. (R2) 

•^0 Jo \=i ^ k=l 

Thus, TZ^f\X) = r['^^(A; 1). Laplace transform {t s) of the above expression yields 

/oo poo 
d\2---j^ rfA^C^|A^({A})inA7e-^"^ (R3) 

This, on using ([2]) and ([5]), gives 

^f\Ms) = ^^<'(^A) = r(^),,'-^<'(.A). (B.4) 

Thus the level density for fixed trace ensemble can be obtained from the level density 
of ordinary Wishart-Laguerre ensemble by taking inverse Laplace transform (s — )■ t) of 
the expression on right hand side of ( 1B.4I) . and by substituting t = 1. 

Appendix C. Proofs of equations (I3ip — (1351) 

We consider the expression 

Rf\N; a;x)= (.) + H'tf {x)L'^.f' {x)^L'i'.f (x)L<;°«' (x)] (C.l) 

for level density given in (CD. Note that Lf\x) = for j < 0. The above expression 
(and other equivalent ones) can be obtained from ((Tj) using the results given in [50] . We 
now use the following standard expansion [50] of Laguerre polynomial in ( IC.ll) : 



^ r(6 + /i + i)r(n -12 + i)r(/x + 1) ' ^ 
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This gives an expression for the level density ([7]) involving double sum. Substituting 
this expression in f lB.4|) and performing the inverse Laplace transform using the result 



gives the result for level density of Schmidt eigenvalues. This result is a double sum, 
however one of the sums can be performed in terms of hypergeometric function defined 
in fl28|) . thus leaving the final result fl3Tl) as a single sum. 

To derive ( 132|) . ( I33l) we use again ( ]C.1|) and the expansion ( ]C.2I) for Laguerre 
polynomial in (|8]) and ([9]). Along with (]C.3|) we also need the following Laplace inverse 
result: 

C-^[s^e'^T{a + 1, sX)] = {t - 2A)-"-^-^A° 2^1 (^1, -a; -a - 7; 2 - e(t - 2A).(C.4) 
Equation f l35|) can similarly be derived from f[T2|) . 

Appendix D. Auxiliary gamma function intergral method 

We introduce an auxiliary gamma function intergral in ( HOl) as 



00 

e ■ r" 



/■oo />oo f \ 



r(si) J, 

- 1 \An{{\WWX"j d\i...d\N dr. (D.l) 
This can be reduced to the following after some manipulations: 



00 

e ' r" 



poo poo f N \ f ^ \ 



X I Ajv({A}) 1^ W A^ d\i...d\N dr. (D.2) 

i=i 

where is the digamma function as defined in fH51) . Now substituting Aj = Xj/r, we 
obtain 



CPOO 
N -r n- 



m) 

Ajm)ff](^y^...^dr 

\ I r I / \ J-J-Vt^/ r r 



■f...f(E?...)<|;?-.) 



.k=l J \i=\ 



'0 r 2 JO 

Af 

X |Aiv({a;})|^ JJx5^da;i...cixivcir. (D.3) 
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The choice Q = f3MN/2 + 1 gives 



TV 

|Ajv({x})|'' JJ^x^ dxi...dxN dr 




j=i 



V 2 -J T(eM!l + 



foo poo / ^ \ 

i)X -I iS"^'""v 



2 

x|Aiv({x})|^ JJx^e"''^ dxi...dxN 



N 



poo poo I ^ \ 

/ ■■■/ V'a;fclnxfc P({a;})(ia;i...cixAr. 



(D.4) 

Using the ratio of constants, Cm/Cn, from (]A.7P we eventually obtain (142|) . 



References 

[1] Horodecki R, Horodecki P, Horodecki M and Horodecki K 2009 Quantum entanglement Rev. Mod. 
Phys. 81 865 

[2] Reid M D, Drummond P D, Bowen W P, Cavalcanti E G, Lam P K, Bachor H A, Andersen U L 
and Leuchs G 2009 The Einstein-Podolsky-Rosen paradox: From concepts to applications Rev. 
Mod. Phys. iSl 1727 

[3] Einstein A, Podolsky B, Rosen N 1935 Can quantum-mechanical description of physical reality be 

considered complete? Phys. Rev. 147 7771 
[4] Schrodinger E 1935 Discussion of probability relations between separated systems Math. Proc. 

Camb. Philos. Sac. 31 555 
[5] Schrodinger E 1936 Probability relations between separated systems Math. Proc. Camb. Philos. 

Soc. \32im 

[6] Bennett C H and Wiesner S J 1992 Communication via one- and two-particle operators on Einstein- 
Podolsky-Rosen states Phys. Rev. Lett. 69 2881 
[7] Harrow A, Hayden P and Leung D 2004 Superdense coding of quantum states Phys. Rev. Lett. 

192 i87gnn 

[8] Bennett C H ei al. 1993 Teleporting an unknown quantum state via dual classical and Einstein- 
Podolsky-Rosen channels Phys. Rev. Lett. 70 1895 

[9] Bennet C H and Brassard G 1984, in Quantum Cryptography: Public Key Distribution and Coin 
Tossing, Proceedings of the IEEE International Conference on Computers, Systems & Signal 
Processing (CS Press), 175 
[10] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge 
University Press). 

[11] Lubkin E 1978 Entropy of an n-system from its correlation with a fc-reservoir J. Math. Phys. 
19 1028 

[12] Page D N 1993 Average entropy of a subsystem Phys. Rev. Lett. 1 71 12"9n 
[13] Sen S 1996 Average entropy of a quantum subsystem Phys. Rev. Lett. 177 II 



19 



[14] Zyczkowski K and Sommers H-J 2001 Induced measures in the space of mixed quantum states J. 

Phys. A: Math. Gen. 34 7lTT 
[15] Cappellini V, Sommers H-J, and Zyczkowski K 2006 Distribution of G concurrence of random pure 

states Phys. Rev. A \74 0623221 
[16] Giraud O 2007 Distribution of bipartite entanglement for random pure states J. Phys. A: Math. 

Theor. 40 2793 

[17] Kubotani H, Adachi S and Toda M 2008 Exact formula of the distribution of Schmidt eigenvalues 
for dynamical formation of entanglement in quantum chaos Phys. Rev. Lett. 1100 2405011 

[18] Adachi S, Toda M and Kubotani H 2009 Random matrix theory of singular values of rectangular 
complex matrices I: Exact formula of one-body distribution function in fixed-trace ensemble 
Ann. Phys. 324 22781 

[19] Vivo P, Majumdar S N, and Bohigas O 2007 Large deviations of the maximum eigenvalue in 

Wishart random matrices J. Phys. A: Math. Theor. 40 4317 
[20] Majumdar S N, Bohigas O, and Lakshminarayan A 2008 Exact minimum eigenvalue distribution 

of an entangled random pure state J. Stat. Phys. 131 33] 
[21] Vivo P 2010 Entangled random pure states with orthogonal symmetry: exact results J. Phys. A: 

Math. Theor. 43 405206 

[22] Majumdar S N 2010 Extreme eigenvalues of Wishart matrices: Application to entangled bipartite 
system larXiv:1005.4515i 

[23] Nadal C, Majumdar S N, and Vergassola M 2010 Phase transitions in the distribution of bipartite 

entanglement of a random pure state Phys. Rev. Lett. 104 110501 
[24] Nadal C, Majumdar S N, and Vergassola M 2011 Statistical distribution of quantum entanglement 

for a random bipartite state J. Stat. Phys. 142 403 
[25] Mehta M L 2004 Random Matrices (New York: Academic) 

[26] Forrester P J 2010 Log-Gases and Random Matrices (Princeton University Press) 

[27] Wishart J 1928 The generalised product moment distribution in samples from a normal 

multivariate population Biometrika 20A 32 
[28] Zyczkowski K and Sommers H-J 2001 Induced measures in the space of mixed quantum states J. 

Phys. A: Math. Gen. m 71111 
[29] Bandyopadhyay J N and Lakshminarayan A 2002 Testing statistical bounds on entanglement using 

quantum chaos Phys. Rev. Lett. 89 060402 
[30] Kubotani H, Toda M, Adachi S 2006 Universality in dynamical formation of entanglement for 

quantum chaos Phys. Rev. A 74 0323141 
[31] Slevin K and Nagao T 1994 Impurity scattering in mesoscopic quantum wires and the Laguerre 

ensemble Phys. Rev. B 50 2380 
[32] Verbaarschot J 1994 Spectrum of the QCD Dirac operator and chiral random matrix theory Phys. 

Rev. Lett. 72 2531 

[33] Ambjorn J, Makecnko Y, and Kristjansen C F 1994 Generalized Pcnner models to all genera Phys. 
Rev. D 50 5193 

[34] Akemann G 1997 Universal correlators for multi-arc complex matrix models Nucl. Phys. B 507 4751 
[35] Bouchaud J-P and Potters M 2001 Theory of Financial Risk and Derivative Pricing (Cambridge: 

Cambridge University Press). 
[36] Burda Z and Jurkiewicz J 2004 Signal and noise in financial correlation matrices Physica A \344 67l 
[37] Telatar I E 1999 Capacity of multi-antenna Gaussian channels Euro. Trans. Telecommun. 110 5851 
[38] Kumar S and Pandey A 2010 Random matrix model for Nakagami-Hoyt fading LEEE Trans. Inf. 

Th. 56 2360 

[39] Alter O, Brown P O, Botstein D 2000 Singular value decomposition for genome- wide expression 

data processing and modeling PNAS 97 10101 
[40] Hoher N S, Mitra M, Maritan A, Cieplak M, Banavar J R and Fedoroff N V 2000 Fundamental 

patterns underlying gene expression profiles: Simplicity from complexity PNAS 9 7 84091 
[41] Edelman A 1989 Eigenvalues and condition numbers of random matrices |P/t.£>. Dissertation, MIT 



20 



[42] Kumar S and Pandey A 2009 Universal spectral correlations in orthogonal-unitary and symplectic- 

unitary crossover ensembles of random matrices Phys. Rev. E 79 026211 
[43] Dumitriu I and Edelman A 2002 Matrix models for beta ensembles /. Math. Phys. 143 58301 
[44] Dyson F J 1972 A class of matrix ensembles J. Math. Phys. 113 901 

[45] Nagao T and Wadati M 1991 Correlation functions of random matrix ensembles related to classical 

orthogonal polynomials J. Phys. Soc. Jpn. 60 3298 
[46] Adler M, Forrester P J, Nagao T and Moerbeke P van 2000 Classical skew orthogonal polynomials 

and random matrices J. Stat. Phys. 199 1411 
[47] Pandey A and Ghosh S 2001 Skew-orthogonal polynomials and universality of energy- level 

correlations Phys. Rev. Lett. 87 024102 
[48] Ghosh S and Pandey A 2002 Skew-orthogonal polynomials and random-matrix ensembles Phys. 

Rev. E \65 0462211 

[49] Kumar S and Pandey A 2011 Crossover ensembles of random matrices and skew-orthogonal 

polynomials Ann. Phys., in press 
[50] Szego G 1975 Orthogonal polynomials (Providence, R I: American Mathematical Society). 
[51] Marcenko V A and Pastur L A 1967 Distribution of eigenvalues for some sets of random matrices 

Math. USSR-Sbormk 1 457 
[52] Vinayak and Pandey A 2010 Correlated Wishart ensembles and chaotic time series Phys. Rev. E 

81 036202 

[53] Zanardi P, Zalka C, and Faoro L 2000 Entangling power of quantum evolutions Phys. Rev. A 
62 030301 (RO 

[54] Scott A J 2004 Multipartite entanglement, quantum-error-correcting codes, and entangling power 

of quantum evolutions Phys. Rev. A 69 052330 
[55] Giraud O and Georgeot B 2005 Intermediate quantum maps for quantum computation Phys. Rev. 

A 72 042312 

[56] Weinstein Y S and Hellberg C S 2005 Entanglement generation of nearly random operators Phys. 
Rev. Lett. 95 030501 

[57] Znidaric M 2007 Optimal two-qubit gate for generation of random bipartite entanglement Phys. 
Rev. A 76 0123181 



